Submitted to ApJ Letters. Oct. 11, 1999 

Preprint typeset using lAT^^X style cniulatcapj v. 04/03/99 



> 

o 



> 

0^ 



COMPTON DRAGGED GAMMA-RAY BURSTS ASSOCIATED WITH SUPERNOVAE 

Davide Lazzati^, Gabriele Ghisellini 

Osservatorio Astronomico di Brcra, Via Bianchi 46, 1-23807 Merate (Lc), Italy 

Annalisa Celotti 

Scuola Internazionale Superiore di Studi Avanzati, via Beirut 2/4, 1-34014 Trieste, Italy 

AND 

Martin J. Rees 

Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, U.K. 
Submitted to ApJ Letters, Oct. 11, 1999 

ABSTRACT 

It is proposed that the gamma-ray photons that characterize the prompt emission of Gamma-Ray 
Bursts are produced through the Compton drag process, caused by the interaction of a relativistic fireball 
with a very dense soft photon bath. If gamma-ray bursts are indeed associated with Supernovae, then the 
exploding star can provide enough soft photons for radiative drag to be effective. This model accounts 
for the basic properties of gamma-ray bursts, i.e. the overall energetics, the peak frequency of the 
spectrum and the fast variability, with an efficiency which can exceed 50%. In this scenario there is no 
need for particle acceleration in relativistic coUisionless shocks. Furthermore, though Poynting flux may 
be important in accelerating the outflow, no magnetic field is required in the gamma-ray production. 
The drag also naturally limits the relativistic expansion of the fireball to F < 10'*. 

Subject headings: gamma rays: bursts — supernovae: general — radiation mechanisms: non-thermal 
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1. INTRODUCTION 

In the leading scenario for GRBs and afterglows, the 
gamma-ray event is produced by internal shocks in a 
hyper-relativistic inhomogeneous wind (Rees & Meszaros 
1994) while the afterglow is produced as the fireball drives 
a shock wave in the external interstellar medium (Meszaros 
& Rees 1997). Even if there is a large consensus that 
both gamma-rays and afterglow photons are produced by 
the synchrotron process, recently some doubts have been 
cast on the synchrotron interpretation for the burst itself 
(Liang 1997; GhiseUini & Celotti 1998; GhiseUini, Celotti 
& Lazzati 1999). 

The nature of the progenitor is still a matter of active 
debate, as the sudden release of a huge amount of energy 
in a compact region generating a fireball does not keep 
trace of the way this energy has been produced. For this 
reason, the study of the interactions of the fireball with 
the surrounding medium seems to be the most powerful 
mean to unveil the GRB progenitor. At least two models 
are in competition: the merging of a binary system com- 
posed of two compact objects (Eichler et al. 1989) and 
the Hypernova-CoUapsar model (Woosley 1993, Paczyhski 
1998), i.e. the core collapse of a very massive star to form 
a black hole. 

After the discovery and the multiwavelength observa- 
tions of many afterglows, circumstantial evidence has ac- 
cumulated for GRB exploding in dense regions, associated 
to supernova-like phenomena. In fact, (a) host galaxies 
have been detected in many cases (Sahu et al. 1997; see 
Wheeler 1999 for a review), and some of them show star- 
burst activity (Djorgovski et al. 1998, Hogg & Fructher 
1999); (b) large hydrogen column densities have sometimes 



been detected in X-ray afterglows (Owens et al. 1998); (c) 
non-detections of several X-ray afterglows in the optical 
band can be due to dust absorption (Paczyhski 1998); (d) 
a possible iron line feature has been detected in the X-ray 
afterglow of GRB 970508 (Piro et al. 1999, Lazzati et al. 
1999) and (e) the rapid decay with time of several after- 
glows can be explained by the presence of a pre-explosion 
wind from a very massive star (Chevalier & Li 1999). More 
recently, the possible presence of supernova (SN) emission 
in the late afterglows light curves of GRB 970228 (Galama 
et al. 1999, Reichert 1999) and GRB 980326 (Bloom et al. 
1999) has added support in favor of the association of some 
GRBs with the final evolutionary stages of massive stars. 

Although in these models the available energy is larger 
than in the case of compact binary mergers, the very small 
efficiency of internal shocks (see, e.g., Spada, Panaitescu & 
Meszaros 1999) seems to be inconsistent with the fact that 
more energy can be released during the burst proper than 
the afterglow (Paczyhski 1999, see also Kumar & Piran 
1999). 

In this letter we show that if GRBs are associated with 
supernovae, Compton drag inside the relativistic wind can 
produce both the expected energetics and the peak en- 
ergy of the spectrum of a classical long duration GRB. In 
this new scenario the efficiency is not limited by internal 
shock interactions, and the successful modeling of after- 
glows with external shocks is left unaffected. 

The Compton drag effect has been already invoked for 
GRBs by Zdziarski et al. (1991) and Shemi (1994). Cos- 
mic background radiation (at high redshift), central re- 
gions of globular clusters and AGNs were identified as 
plausible sources of soft photons, but none of these scenarii 
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was able to account for all the main properties of GRBs. 
However, the growing evidence of association of GRB ex- 
plosions with star-forming regions and supernovae opens 
new perspectives for this scenario. 

2. COMPTON DRAG IN A RELATIVISTIC WIND 

We consider a relativistic (F 3> 1) wind of plasma prop- 
agating in a bath of photons with typical energy Csood- A 
fraction ^ min(l,TT) of the photons are scattered by the 
inverse Compton (IC) effect to energies e ^ F^esocd, where 
tt is the Thomson opacity of the wind. Due to relativistic 
aberration, the scattered photons propagate in a narrow 
cone forming an angle l/F with the velocity vector of wind 
propagation. By this process, a net amount of energy Ecb 
is converted from kinetic energy of the wind to a radiation 
field propagating in the direction of the wind itself, where 
Ecu ~ niin(l,rT) y Wrad (r^ - !)■ V is volume filled by 
the soft photon field of energy density Uiad swept up by 
the wind. 

Let us assume that the GRB fireball, instead of being 
made by a number of individual shells (see e.g. Lazzati et 
al. 1999), is an unsteady (both in velocity and density) 
relativistic wind, expanding from a central point. After 
an initial acceleration phase, the density of the outflowing 
wind decreases with radius as n(r) oc r~^, giving a scatter- 
ing probability ~ min[l, (r/ro)~^], where tq is the radius 
at which the scattering probability equals unity. After the 
first scattering, the photons propagates in the same direc- 
tion of the flow and the probability of a second scattering 
is reduced by a factor ^ F^. 

If such a wind flows in a radiation field with energy den- 
sity Mrad(?'), the total energy transferred to the photons 
when the fireball reaches a distance R is given by^: 
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where for simplicity we assume that a constant F has been 
reached (see also Section 3). The transparency radius rg 
depends on the baryon loading of the fireball, which is pa- 
rameterized by 77b = E/{Mc^), where E/M is the ratio 
between the total energy and the rest mass of the fireball. 
Then tq is given by"^: 



ro = 5.9 X 10"£:,\/%7u^(? cm. 



2.1. 
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A simple scenario 
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We initially consider a simple scenario which can il- 
lustrate the basic features of the Compton drag effect. 
Let us assume that the GRB is triggered at a time Ai 
(of the order of a few hours) after the explosion of a 
supernova (Woosley et al. 1999; Cheng & Dai 1999). 
By this time, the supernova ejecta, moving with veloc- 
ity /3snc, have reached a distance i?sN = I'SNAi ~ 5.4 x 
10"'^'^/3{gN,-i}(Ai/5 hr) cm. Let us also imagine that the 
supernova explosion is asymmetric, e.g. with no ejecta in 
the polar directions. Despite this asymmetry, the ejecta 



uniformly fill with radiation the entire volume within i?sN • 
If i?sN < fQ, the energy extracted by Compton drag is: 
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i?SN > ro. (4) 



'^, with a black body emission at 
10^ K. It follows that in this case 



According to Woosley et al. (1994), the average lumi 

nosity of a type II supernova'' during At is of the order 

of LsN ^ 10^"* erg s 

a temperature Tsn ^ 

Urad == aT^ 

i?SN assumed above 

in extracting the fireball energy is very large; from Eq, 

we obtain: 
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Note here that a high efficiency can be reached even for 
F ~ 100. Note that the drag itself can limit the maximum 
speed of the expansion - even in a wind with a very small 
barion loading - as discussed in Sect. & Each seed pho- 
ton is boosted by ~ 2F^ in frequency, yielding a spectrum 
peaking at hv 
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2F^(3fcrsN) - 0.5FiT{SN.5} MeV. 
. A more realistic scenario 



The previous scenario requires that the GRB explodes a 
few hours after a supernova. There is however a plausible 
alternative, independent of whether the massive (> 3OM0) 
star (assumed to be the progenitor of the GRB) ends up 
with a supernova explosion or not, and can produce a 
gamma-ray burst even if the relativistic flow and the core 
collapse of the progenitor star are simultaneous or sepa- 
rated by a relatively small time interval (Woosley et al. 
1999; MacFadyen, Woosley & Heger 1999). 

In fact there is a somewhat general consensus (e.g. Mac- 
Fadyen & Woosley 1999; Aloy et a. 1999, but see also 
Khokhlov et al. 1999) that a relativistic wind can flow 
in a relatively baryon free funnel created by a bow shock 
following the collapse of the iron core of the star. Even 
if details of this class of models are still controversial, the 
formation of the funnel seems to be a general outcome. Let 
us estimate its luminosity and more precisely the amount 
of energy in radiation crossing the funnel walls at a time 
if after its creation. With respect to the total luminos- 
ity of the star, assuming it radiates at its Eddington limit 
~ LEdd, there would be a reduction by the geometrical fac- 
tor equal to the ratio of the funnel to star surfaces, which 
is of the order of the funnel opening angle -d. However, im- 
mediately after its creation, the funnel luminosity is much 
larger than "iJLEdd, due to two effects which we discuss in 
turn. 

First the walls of the funnel contain an enhanced amount 
of radiation with respect to the surface layers of the star: 
the radiation once "trapped" in the interior of the star 
can escape through the funnel walls, thus enhancing the 
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luminosity inside the funnel for a short time. Photons 
produced at a distance s from the wall surface cross it at 
a time t{ ~ Tss/c = ans'^ /c, where a is the relevant cross 
section. This compares with the Kelvin time ix ~ anR^/c 
needed for radiation to reach the star surface, yielding 
s/Ri, ^ (tf/tK)^/^. After the time if, the radiation pro- 
duced in the layer of width ds crosses the funnel surface 
carrying the energy dEf '^ •dTi.L-Eddds/c, the correspond- 
ing luminosity being: 
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For if = 100 s and a IOMq star with i?* ~ 10^" cm 
(r^ w 10^), this effect can enhance the funnel luminosity 
by ~ 10^. 

Let us now consider a second plausible enhancing fac- 
tor. If the funnel has been produced by the propagation 
of a bow-shock in the star, the matter in front of the ad- 
vancing front is compressed, with a pressure increase of 
7W^, where A4 is the Mach number of the shock in the 
star. This (optically thick) gas then flows along the sides 
of the funnel and relaxes adiabatically to the pressure of 
the external matter (its original pressure). The result is 
that the funnel is surrounded by a sheath (cocoon) with 
density lower than that of the unshocked stellar material 
by a factor TW^^^ (a polytropic index of 4/3 has been used 
in the adiabatic cooling) . The diffusion of photons through 
this rarefied gas into the funnel is then even faster, result- 
ing in a further increase of the luminosity by A4^^'^ ~ 200, 
where a shock speed /3sfC — 0.1c (MacFadyen & Woosley 
1999) and a sound speed /3sC = 10~''c have been assumed. 

By taking into account both effects, the funnel luminos- 
ity corresponds to: 
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which leads to an energy loss for Compton drag Lqd ^ 
F^Lf - lO''9t9_iF|(M*/lOM0), to be compared with the 
observed luminosity (Lgrb) "^ lO^^Tn??.]^ erg s~^. Here 
the average luminosity is considered over the entire burst 
duration: for single pulses, we should take into account an 
extra factor F^ in Eq. H due to the Doppler contraction of 
the observed time. 

The typical radiation temperature associated with this 
luminosity, assuming a black body spectrum, is enhanced 
with respect to the temperature of the star surface by 
[L(/{^LEdd)V^^ r- (r.i?.)i/«(cif)-i/8(;j^j/^^)3/i6. Adopt- 
ing the numerical values used above, the enhancement is 
of the order of 50, corresponding to a funnel temperature 
Tf ^ 2 X 10^ K (for a surface temperature of the star of 
^5000 K). This value is similar to the one estimated in the 
simple scenario of the previous subsection and thus leads 
to similar Compton frequencies. 

3. PROPERTIES OF THE OBSERVED BURSTS 

If the wind is homogeneous the spectrum of the scat- 
tered photons resembles that of the incident photons, i.e. 
a broad black-body continuum peaked at a temperature 
Tdrag ~ 2F^T. While the observed characteristic photon 
energy would be therefore e ~ 0.5F2T5(1 -I- z)~^ MeV, 



in good agreement with the observed distribution of peak 
energies of BATSE GRBs (assuming again F — 100, see 
below), the spectrum would not reproduce the observed 
smoothly broken power-law shape (Band et al. 1993). 
The assumptions of a perfectly homogeneous wind and of 
an isothermal radiation field are however very crude, and 
one might reasonably expect that different regions of the 
wind are characterized by different values of F and dif- 
ferent soft field temperatures. If we assume, e.g., that 
the temperature of the soft photon field varies with ra- 
dius according to a power-law T{r) oc r~ , the time in- 
tegrated spectrum will have a high energy power-law tail 
F{h') oc v 5~ . In addition, the bulk Lorcntz factor of 
the flow is likely to vary on a timescale much shorter than 
the integration time required to obtain a spectrum with 
the BATSE data (~ 1 s), and hence the analysed spec- 
tra are the superposition of drag spectra by many differ- 
ent Lorentz factors. A third effect adding power to the 
high energy tails of the spectrum is the reflection of up- 
scattered photons in the pre-supernova wind. This pho- 
tons are scattered again by the fireball and can reach ener- 
gies of ~ 0.5F MeV - SOFz MeV. The computation of the 
actual spectrum resulting from all these effects depends 
from many assumptions and is beyond the scope of this 
work. 

The effects described above, which can increase the fun- 
nel luminosity over the Eddington limit, take place in non- 
stationary conditions. At the wind onset, it is likely that 
the temperature gradient in the walls of the funnel is large, 
but this is soon erased due to the high luminosity of the 
walls. This causes both the total flux and the character- 
istic frequency of the soft photons to decrease, and hence 
a hard-to-soft trend is expected. Moreover, it has been 
shown by Liang & Kargatis (1996) that the peak frequency 
of the spectrum in a single pulse at time i is strongly re- 
lated to the flux of the pulse integrated from the beginning 
of the pulse to the time i. In our scenario, this behaviour 
can be easily accounted for if we consider a shell slowed 
down by the drag itself: the Lorentz factor (and hence the 
peak frequency of the spectrum) at a time i is related to 
the energy lost by the shell, i.e. to the integral of the flux 
from the beginning of the pulse to the time i. 

The observed minimum variability time-scale is related 
to the typical size of the region containing the dense seed 
photon field, which corresponds to either Ri, or i?sN de- 
pending on which of the two scenarios described above 
applies. The relevant light crossing time - divided by the 
time compression factor - is thus 
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Longer time-scales are instead expected if the relativistic 
wind is smooth and continuous. 

Another interesting feature of this scenario is the pos- 
sibility that the bulk Lorentz factor of the wind is self- 
consistently limited by the drag itself. The pressure of 
the soft photons starts braking the fireball in competi- 
tion with the pressure of internal photons. The limiting 
Lorentz factor is hence reached when the internal pressure 
p'fYj oc (To/r)"* is balanced by the pressure of the exter- 
nal photons as observed in the fireball comoving frame 
p' oc F^ Tgj^ (1 -I- tt)"^, where tt is the scattering optical 



depth of the wind. This gives: 
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where i?o is the radius at which the firebaU is released. 
Equation | reduces to Tum -^ 10''(To,ii/rsN,5)^/^ if the 
firebaU becomes transparent before reaching the coasting 
phase. With such high T the Compton drag would be max- 
imally efficient , causing the fireball to immediately decel- 
erate until its r reaches the value given by Lcd = LfV^ , 
implying: 
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These limits are in general smaller than the maximum F 
set by the baryon load only, but still in agreement with 
values recently inferred for GRB 990123 (Sari & Piran 
1999). In addition, it is likely that the external parts of 
the relativistic wind, which are in closer connection with 
the funnel walls, are dragged more efficiently then the cen- 
tral ones, since at the beginning the soft photons coming 
from the walls can penetrate only a small fraction of the 
funnel before being up-scattered by relativistic electrons. 
This may result in a polar structured wind, with higher 
Lorentz factors along the symmetry axis, gradually de- 
creasing as the polar angle increases. 

4. DISCUSSION 

A crucial requirement of our model is the association of 
GRBs with the final evolutionary stages of very massive 
stars, as these provide the large amount of seed photons 
emitted at distances ~ 10^3 cm from the central trigger, 
which are neeeded for the Compton drag to be efficient. 

The efficiency of conversion of bulk kinetic energy of 
the flow into gamma-ray photons is large, solving the ob- 
servational challenge of gamma-ray emission being more 
energetic than the afterglow one (Paczyhski 1999). Fur- 
thermore in this scenario there is no requirement for ei- 
ther efficient acceleration in collisionless shocks or the 



presence/generation of an intense (equipartition) magnetic 
field, although Poynting flux may still be important in ac- 
celerating the outflow (being more efficient than neutrino 
reconversion into pairs). 

We have investigated the main properties of a GRB 
produced by Compton drag in a relativistic wind in a 
very general case. A moderately beamed burst [d <^ 10°, 
Woosley et al. 1999) can be thus produced and, without 
any fine tuning of the parameters, the basic features of 
classic GRBs are accounted for. 

In particular, the peak energy of the burst emission sim- 
ply reflects the temperature of the supernova seed photons, 
up-scattered by the square of the bulk Lorentz factor. 
The simplest hypothesis predicts a quasi-thermal spec- 
trum, however it is easy to imagine an effective multi- 
temperature distribution which would depend on uncon- 
strained quantities such as the variation of the spectrum 
of the SN photons with radius and the degree of inhomo- 
geneity of the wind. 

Although in this scenario there is no requirement for 
internal shocks to set up, they can of course occur, con- 
tributing a small fraction of the observed gamma-ray fiux. 
On the other hand, the wind is expected to escape from 
the funnel of the star with still highly relativistic motion, 
so that an external shock can be driven in the interstellar 
medium and produce an afterglow, similar to the scenario 
already studied by several authors. It is likely that this af- 
terglow would develop in a non-uniform density medium, 
due to the presence of the massive star wind occurring 
before the supernova explosion (Chevalier & Li 1999). 
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